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Abstract 

We discuss the global symmetries of the high temperature phase of QCD with Nj 
massless quarks. We show that the U{Nf) x U{Nf) symmetries are only violated by 
operators of dimension > 3Nf. For Nf > 2 this implies that the thermal two-point 
correlation functions of the rj' and vr'^'s are identical. We discuss the implications of 
this for the chiral phase transition at finite temperature. 
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Recently, Cohen has used QCD inequalties to argue that the high temperature phase 
of massless QCD with Nj flavors of quarks is effectively symmetric under a global U{Nf) x 
U{Nf) symmetry rather than just SU{Nf) x SU{Nf). Cohen examines the two-point cor- 
relation functions of different operators such as the rj' {q'j5q) and the vr" {qT""y5q) which are 
in the same U{Nf) x U{Nf) multiplet but not in the same SU{Nf) x SU{Nf) multiplet. He 
argues that in the massless (or chiral) limit: —>■ 0, the difference between the respective 
two-point correlation functions approaches zero. 

In this letter we show that the U{Nf) x U{Nf) symmetry is not completely restored in 
the high temperature phase although its breaking can only be manifested in operators of 
dimension > 3A^/. Thus the U{Nf) x U{Nf) symmetry is only restored for the two-point 
correlators when Nf > 2. We concentrate on the relationship between (spontaneous) chiral 
symmetry breaking, the axial anomaly and topological charge and clarify some subtle points 
in the argument of [|1|. 

It is well-known that the solution of the U{1)a problem requires 

(A) the axial anomaly relation (m^ = 0) 

d,4 = #tr[FF], (1) 

and 

(B) the presence of gauge configurations with non-zero topological charge 

j d'xtr[FF]. (2) 

We can understand this by examining the behavior of correlators under U{1)a transfor- 
mations. Let| 

{O) = ^ I D[A] e"^^*^ J D^Di^e-I'f'^'^'^''^^ O (3) 

where O is an operator built out of ip^tl) and gauge fields. Now consider O after an axial 
transformation: 

O^O^ = O (V',V'') (4) 

where 

^ ^ ^' = e^°T57/> (5) 

^ ^ ^' = ^e^"^^ (6) 

^We assume throughout this paper a non-perturbative regulator which conserves chiral symmetry and is 
controlled by a UV scale A. See, e.g., |^. 
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The change in the expectation of O is given by 

6^{o) = (a) - {o) 

= ^ j D[A]e-^''^' J Dt/jDil^e-I^^'^'^''^^ [Oa-O]. (7) 

The integral in (|7]) can be evaluated by a change of variables: ip ^ ip', ip ^ ip'. The only 
subtlety is that the Jacobian for the change of variables induces the anomaly factor in the 
measure of the functional integral IQ. In the chiral limit, this yields 

6^{0) = ( [e^-^ -1]0) = i [e'^"' - 1] (§) 

where ( )i, denotes the expectation value taken in the sector of the gauge field configuration 
space with topological charge u. Thus we see that the physical effect of the U{1)a anomaly 
is only manifested in sectors of the functional integral with non-zero topological charge. 

The operator relation in (0), which we used to derive (||) follows from the ultraviolet 
(UV) behavior of the theory and is therefore unaffected by temperature. One can easily see 
this by repeating the Fujikawa derivation with boundary conditions in the Euclidean time 
direction which are appropriate for T ^ 0. As long as the UV regularization scale is kept 
large compared to T the same result for the anomalous variation of the functional measure is 
obtained. However, we will demonstrate below that above the temperature at which chiral 
symmetry is restored, the contributions of gauge configurations with u ^ to correlators 
of quark operators of dimension < 3A^/ are suppressed (they are in fact a set of 'measure 
zero' in the functional measure) and hence the axial anomaly has no effect on the U{1)a 
Ward identities for these operators. For these operators we will show that the right hand 
side of is zero when there is no spontaneous chiral symmetry breaking (i.e. in the high 
temperature phase of QCD) and the current quark masses rriq are taken to zero. 

The presence of massless quarks is known to suppress topological fluctuations. The 
partition function for QCD is 

Z = J D[A] e'^^" Det[p-mg] e'^I^^ 

^ jWA]e^'^'^ (9) 

where Det[Ip— m^] = Hn (^-^n — ^^^q) and the A„ are eigenvalues of the Euclidean Dirac 
equation: pipn = iXn'ipn- We can break Z into contributions from sectors of different winding 
number 

l'=+OC 

z = e''r (10) 

u=—ao 

An index theorem [§] tells us that there must exist a minimum number of zero mode solutions 
when u ^ 0: u = n+ — n^, where n± is the number of right (left)-handed solutions. This 
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implies that in the z/ 7^ sectors Zy must vanish at least as fast as mq ^ in the chiral limit. 
(In this paper we will assume exact SU{N)v isospin symmetry and take the quark mass 
matrix to be proportional to the identity matrix: M = rrigl.) This seems to imply that only 
the i> = sector contributes to Z when massless quarks are present. 

When chiral symmetry breaking is involved, this argument is too naive. This is fortunate 
since we believe that in QCD, which exhibits chiral symmetry breaking, the U{1)a problem 
is indeed solved by the combination of (A) and (B). Heuristically, one might guess that the 
'dynamical mass' acquired by the quarks plays some role in determining whether topological 
fluctuations are suppressed. This issue has been addressed systematically by Leutwyler and 
Smilga 1^ . They find that fluctuations in topological charge are controlled by the parameter 
X = E Vniq, where 



is the chiral symmetry breaking order parameter, and V is the volume of the system. At 
large X topological fluctuations are allowed, whereas as X — > they are suppressed. The 
order of limits in taking V ^ 00 and rriq — is important to the discussion here. If one takes 
rriq — > with V fixed one does not recover the phase of QCD in which chiral symmetries 
are spontaneously broken. This is because symmetries cannot break spontaneously at finite 
volume, and hence for fixed V a non-zero nig is required to bias the vacuum energy and keep 
the system in the broken vacuum. It is only at X >> 1 that finite- volume effects are small 
and we recover the low-energy vacuum of QCD. However in this limit topological fluctuations 
are unsuppressed even though rriq goes to zero. This is essentially a consequence of chiral 
symmetry breaking. 

It is believed that in the high temperature phase of QCD chiral symmetry is restored 
by thermal effects. In other words, at sufficiently high temperature, the minimum of the 
free energy is at S = 0. We can study the theory at finite temperature by imposing (anti-) 
periodic boundary conditions in Euclidean time on the (fermion) boson fields appearing in 
the functional integral. If the period (3 = 1/T is taken to be sufficiently small we will recover 
the high temperature phase. In this phase the subtlety associated with the order of limits 
rriq — > 0, y ^ 00 is no longer importantP] as we are not trying to recover a phase with 
spontaneous symmetry breaking. Instead, we can take rrig before taking the V 00 
(or the UV cutoff A —>■ 00). This has important consequences, as it allows the naive scaling 
arguments in ruq to be used. It also allows us to make our arguments rigorous, as all 
quantities are finite while V and A are kept finite. For example, in the chirally restored 
phase Ziy does indeed vanish like m}q^^^ in the chiral limit. In the absence of dynamical 

^For the order of limits to be important, some physical quantities would have to depend on parameters 
such as, e.g., m^L, where L is the size of our box. However, this is highly implausible in the high-T 
(disordered) phase as there are neither long-range order nor long-range correlations in the heat bath. 




(11) 
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chiral symmetry breaking massless quarks suppress the contributions to Z from sectors of 
non-trivial topology. 

We next consider the contributions from sectors with 7^ to arbitrary correlators in 
the high-temperature phase. By performing the fermionic part of the functional integral, the 
right hand side of (H) can be written as a sum over permutations of functional integrals over 
the gauge field measure with integrands that are products of traces of propagators, gamma 
matrices and flavor matrices (times possibly some functions of the gauge field, which we 
suppress) : 

^a(0) = E E / V [e'"' - 1] n [sA(x„x,)ri5l(xfc,x;)r2...] , (12) 

v^Q perms 

where F is any combination of gamma and flavor matrices. 

In each v sector the measure, /[^/xa] v, approaches zero as where no is the number 
of zero modes and hence no > \v\Nj. The propagators may be written in terms of their 
spectral decomposition 

(13) 

where it is important to remember that Fermi statistics (imposed by the integration over 
Grassmanian variables) forbids any two propagators in the product of sums from sharing an 
eigenvalue Afc. Thus the integrands in (O) diverge at most as (1/mg)" where the operator O 



is of dimension 3n but where the power of divergence must always be < ng. The functional 
integral will never diverge as a result of the — limit since the propagators in O may 
at most 'soak up' all no zero modes []. However, we observe that correlators of all operators 
of dimension < 3A^/ receive no contributions from gauge configurations with i/ 7^ in the 
ruq — >■ limit. The operators whose expectation values in z/ 7^ sectors are non- vanishing 
when niq are precisely the so-called 'tHooft operators ^ induced by instanton processes 
at weak coupling. 

The two-point correlators are of special interest since they determine the number of 
(nearly) massless modes present at and above the high temperature phase transition, which 
is believed to be either second or weakly first order. For example, consider the two-point 
functions for the vr and Vj at finite temperature: 

(r/'(x)r7'(0)) = ( V>,75V'. (x) ^-.Ts^i (0)) (14) 
(7r'^(a:)7r'^(0)) = ( ^r^^^^ [x) ^r^^^^ (0) ) (15) 



"This is a loophole in the argument of |^ where it is argued that if the contribution to an an n-point 
correlator from a particular v sector is non-zero the (n-|-2)-point quark correlator must diverge. If the n- 
point correlator is non-zero all zero modes in that sector have been soaked up by the propagators and the 
introduction of further propagators will not generate any further inverse powers of rtiq. This loophole is 
problematic for the results of |l[ when TV/ = 2. 
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Writing these correlators in terms of exact quark propagators SA{x,y), one finds two types 
of contributions: a disconnected contribution 

i J dfiA tr[rSA{x,x)]Tr[rSA{0,0)] (16) 

and a connected part 

I J dfiA tr[SA{x,0)rSA{x,0)r]. (17) 

Here F = 75 for the 1]' and F = r"75 for the vr". The connected parts ([TTD are identical since 
[t°', Sa] = 0. For the pion, the disconnected part is zero since tr[T°'] = 0. Any r/'-vr'' splitting 



is the result of (ITO) for the rj'. 



The measure in a given u sector goes to zero at least as fast as mJq^^^ as rrig whilst 
the two propagators may soak up only two of the zero modes. Thus for Nf > 2 the difference 
between the correlators vanishes as m ^ and the rj' is degenerate with the massless pions. 
A similar result holds for the entire U{Nf) x U{Nf) multiplet which consists of the 11,6, a 
and 7]' resonances. 

When Nf < 2, however, there are potentially contributions to the correlator difference 
from the z/ = ±1 sector. There is good reason to believe that these contributions are non-zero 
since one-instanton effects which contribute to ([l6|) were found in the weak coupling regime by 
't Hooft . His calculations in the weak coupling approximation are relevant at temperatures 
T ^ Aq(7£), where the effective coupling constant is small and non-perturbative effects can 
be studied in the semiclassical approximation (see also 0). When Nf = 2 instanton effects 
contribute directly to ([T6|), leading to an rj'-n"- mass splitting. It is possible in principle that 
near the chiral phase transition T ~ Tc, where the dilute instanton analysis is not completely 
reliable, that some other effects (such as I-I pair formation [§]) lead to the suppression of 
this splitting, but it is implausible that (|T6|) vanishes completely. 

In the case of Nf = 1 there are potentially contributions to the chiral condensate 



E = y dfiA tr[SAix,x)] 



from the = ±1 sector. These can be seen to be non-zero even at high-temperature 0, and 
combining this with the low temperature analysis of Leutwyler and Smilga |Q , one reaches 
the conclusion that the U{1)a violating chiral condensate stays non-zero for all temperatures 
with U{l)v being the only unbroken global symmetry. 

We end with a summary of our conclusions for different values of Nf, and the correspond- 
ing implications for the chiral phase transition. 

• Nf = 1: There is no chiral phase transition, (ipip) remains non-zero at high tem- 
peratures, and the dynamically generated quark mass decreases smoothly to zero as 
T 00. 
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Nf = 2: The SU (2) x SU (2) global symmetry is restored in the high temperature phase. 
The tj'-tt"' splitting is non-zero, but decreases smoothly to zero with temperature as 
determined by the instanton density: m^/^ ~ (^)'^- At T ~ Tc the r]' is massive and the 
symmetry remains SU{2) x SU{2). Renormalization group calculations ^ (based on 
the e-expansion) in an effective linear sigma model with this symmetry indicate that 



the phase transition is second order and QCD lattice simulations |]T0[ appear to be in 
agreement with this prediction. It is worth mentioning that if the rj'-ir"' splitting is 
non-zero but small at T ~ Tc, the effective model for the system as T ^ Tc from above 
may have an approximate U{2) x U{2) symmetry. Studies of U{2) x U{2) dynamics 
P generically lead to a fluctuation- induced first order transition, so it is possible that 
the chiral phase transition occurs through this instability before reaching the second 
order fixed point at length scales much larger than the rj' correlation length. For an 
overview of these and related issues, see 

Nf > 3 (assuming asymptotic freedom): The U{Nf) x U{Nf) global symmetry is 
effectively restored (up to high- dimension operators which are probably irrelevant in 
the IR limit) in the high temperature phase. The 1]' becomes degenerate with the vr^'s 
for T > Tc. The effective models of this chiral phase transition should incorporate a 
U{Nf) X U{Nf) global symmetry. The analysis of the appropriate linear sigma model 
without determinantal interactions [Tl|] suggests that the phase transitions is first 
order with its strength increasing with Nf. 

Finally, it is amusing to note that in QCD theories with Nf just below yiVc, such that 
the theory has a perturbative infra-red fixed point, there is no dynamical chiral sym- 
metry breaking and hence our results also apply there independent of the temperature. 
In such theories the U{Nf) x U{Nf) symmetry applies to all correlators up to some 
high dimension, where there are some exponentially small U{1)a violating interactions 
in the form of high-dimension instanton operators ('tHooft vertices). Strictly speaking 
there are no mesons here since we are at weak coupling. 
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